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An action principle of singular hypersurfaces in general relativity and scalar-tensor type theories 
of gravity in the Einstein frame is presented without assuming any symmetry. The action principle 
is manifestly doubly covariant in the sense that coordinate systems on and off a hypersurface are 
disentangled and can be independently specified. It is shown that, including variation of the metric, 
the position of the hypersurface and matter fields, the variational principle gives the correct set 
of equations of motion: the Einstein equation off the hypersurface, Israel's junction condition in 
a doubly covariant form and equations of motion of matter fields including the scalar fields. The 
position of the hypersurface measured from one side of the hypersurface and that measured from 
another side can be independently variated as required by the double covariance. 



I. INTRODUCTION 

Spatially extended objects such as branes, membranes, shells and walls have been playing many important roles 
in recent progress in theoretical physics including string theory [jy, particle phenomenology |2|-Q, theory of black 
holes cosmology [|]]|] an d so on. Hence, it is important to investigate dynamics of such extended objects. In 

particular, the so called brane world scenario is based on the idea that our four-dimensional universe may be a world 
volume of a brane in a higher dimensional spacetime p[-p|,^0[ . Thus, in the brane world scenario the dynamics of the 
brane is the dynamics of our universe itself and is of the most physical importance. 

It is well-known and is the most commonly adopted picture that the dynamics of extended objects is elegantly 
described as geometrical imbedding of world-volume surfaces into spacetime in a certain limit. In particular, in the 
case of codimension 1, or when the world- volume surface is a hypersurface, the geometrical description becomes 
simpler than other cases with higher codimension. Actually, in general relativity or other theories of gravity in the 



Einstein frame, the classical dynamics of a hypersurface is perfectly described by Israel's junction condition |11|. 

One of the main advantages of the junction condition is that it is manifestly doubly covariant in the sense that 
coordinate systems on and off a hypersurface are disentangled and can be independently specified. More precisely, 
there are three independent coordinate systems: that on the hypersurface, those in two regions separated by the 
hypersurface. In the brane world point of view, the double covariance is important since it allows us to separate the 
coordinate system in our world from that in the higher dimensional spacetime. 

Once the classical dynamics is understood, one would usually like to understand quantum mechanical dynam- 
ics |T^-|2l|]. For this purpose, we would like to obtain the action principle for the system including a hypersurface. 

The easiest way to obtain the action may be to adopt the Gaussian normal coordinate system based on the 
hypersurface and to consider the Einstein-Hilbert action with a delta function source. The action obtained in this 
way gives the correct set of equations in the coordinate system, provided that the position of the hypersurface and 
coordinates in a neighborhood of the hypersurface are fixed by the Gaussian normal coordinate condition. However, 
in this method we lose the double covariance: coordinates on the hypersurface is a part of coordinates off the 
hypersurface so that the coordinates satisfy the Gaussian normal coordinate condition. The loss of the double 
covariance is regrettable. 

Actually, as far as the author knows, a doubly covariant action principle has not yet been obtained in the literature. 
One of the main difficulties seems due to the fact that the spacetime metric on one side of the hypersurface and 
that on another side are independent variables in the variational principle. Hence, a question arises: How can we 
ensure the regularity of the intrinsic geometry of the hypersurface without entangling the coordinate systems on and 
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off the hypersurface? This question will be answered in this paper as a manifestly covariant action principle will be 
presented. 

Another difficulty is due to the fact that the double covariance requires inclusion of the position of the hypersurface 
as a dynamical variable in the action principle. Actually, in the doubly covariant formulation of the junction condition, 
it is easy to see that variables specifying the position are not invariant under coordinate transformation and should 
not be fixed J22|. More about why we need to include the position of the hypersurface will be explained from the 
brane world point of view in Sec. M. Here, we mention that, since coordinate systems in two sides of the hypersurface 
are independent, the position of the hypersurface measured from one side and that measured from another side should 
be independently variated in the variational principle. 

It may be worth while reviewing the present status in the literature regarding the second difficulty. However, 
the author knows only a few papers referring this point. Here, we only quote a sentence from one of them: ' The 
variational equations that arose from the unreduced Hamiltonian action were not strictly consistent in a distributional 
sense, but we were able to localize the ambiguity into the single equation that arises by varying the action with respect 
to the shell position ^0|'. One might think that another paper had obtained the correct set of equations, but in 
that paper the position of the hypersurface measured from one side and that measured from another side can not be 
variated independently. Actually, if we would simply variate the position of the hypersurface measured from one side 
and that measured from another side independently, then the variational principle presented in ref. would give 
wrong equations. Moreover, in both of these papers, the hypersurface represents only a dust shell and the coordinate 
systems on and off the hypersurface are not independent. One of them [2^| assumes spherical symmetry, too. 

The purpose of this paper is to provide a manifestly doubly covariant action principle of singular hypersurfaces in 
general relativity and scalar-tensor type theories of gravity in the Einstein frame without assuming any symmetry. 
Besides the scalar fields included in the scalar-tensor type theories, any kind of matter Lagrangian density on the 
hypersurface, which may depend also on the pullback of the scalar fields, can be included. It is shown that, including 
variation of the metric, the position of the hypersurface and matter fields, the variational principle gives the correct 
set of equations of motion: the Einstein equation off the hypersurface, Israel's junction condition in a doubly covariant 
form and equations of motion of matter fields including the scalar fields. As required by the double covariance, the 
position of the hypersurface measured from one side of the hypersurface and that measured from another side can be 
independently variated. 

This paper is organized as follows. In Sec. O a doubly covariant action of a singular hypersurface is derived from 
the standard Einstein-Hubert action. In Secrlinl the variation of the action is calculated for the variations of the 



metric and the position of the hypersurface, and the corresponding equations are obtained. In Sec. IV the variation 
of the action corresponding to the variations of scalar fields are evaluated. Sec. [v| is devoted to a summary of this 
paper and some discussions. 



II. ACTION OF SINGULAR HYPERSURFACE 



Let us consider a D-dimensional spacetime (A4, gMN) and a timelike or spacelike hypersurface S which separates 
A4 into two regions, A4 + and A4-. Since we would like to consider £ as a physical object (eg. the world-volume of 
a brane or the world-volume of a bubble wall in a first-order phase transition) or a physical event (eg. instantaneous 
global phase transition p3f ), we assume that the (D — l)-dimensional intrinsic geometry on £ is regular. On the other 
hand, the _D-dimensional geometry is not necessarily regular on S. 

In the following arguments we shall estimate the action for the system including the singular hypersurface S. We 
assume that the system is described by the action [] 

Itot = I EH + Imatten (1) 

where Ieh is the Einstein-Hubert action with a cosmological constant 

Ieh = ^ f d D x^\g~\{R- 2A), (2) 
and Imatter is the matter action of the form 



For simplicity we do not consider the boundary of M, but it is easy to take it into account by imposing suitable boundary 
conditions and introducing boundary terms appropriate for the boundary condition. 
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Here, {s^f } are D-dimensional coordinate systems in M±, respectively, and is a (D — l)-dimensional coordinate 
system in E. The three coordinate systems can be independent from each other. 

In order to evaluate the gravitational part of the action, we first regularize the Z?-dimensional geometry in a 
neighborhood of E by introducing the finite thickness 5 of the object corresponding to E. Of course, in the final 
step below, we shall take the limit S — > +0, where the hypersurface becomes singular again. Namely, we consider the 
prescription 



Ieh= lim (7| 
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where Mq is a spacetime neighborhood of E representing the regularized object, M± are the two regions separated 
by Mq so that 



Mi D E, M± C M±, lim M 5 ± 

s^+o 



M- 



(5) 



and £>± is the boundary between Mq and -M!j_, respectively. Note that surface terms have been included in Iq ± for 
later convenience but that these exactly cancel each other on common boundaries B±. Each surface term is defined as 
an integral over the (D — l)-dimensional intrinsic coordinates y M on B±, q is the determinant of the induced metric, 
K is the trace of the extrinsic curvature associated with the unit normal n M directed from Mq to M 5 + or from M 5 _ 



to Mq, and e 



gMNn M n N 



±1. 



Next, in order to estimate Iq, we foliate Mq by such a one-parameter family of hypersurfaces E r that So coincides 
with E and that E±i coincides with the boundary B±, respectively. Hence, we can decompose the D-dimcnsional 
Rlcci scalar R as 



R = R (D - 1} + eK 2 - eK^K^ - 2e(Kn M - n M N n N ), M , 



(6) 



where R^ D ^ is the Ricci scalar of the (D — l)-dimensional induced metric on E r , the semicolon represents the 
covariant derivative compatible with gMN, n M is the unit normal to E r directed towards B+, e — gMNn M n N = ±1, 
K^ v is the extrinsic curvature associated with n M , the indices {//, v\ are raised by the inverse of the induced metric, 
and K = Kjf. By integrating over Mq and taking the limit 6 — > +0, we obtain 

(7) 

f 0. We have 



= A / d D x^\g~\ (r {d -^ + eK 2 - eK^K^ - 2A) -> (5 -> +0). 
2k /m* v y 



Here, we have used the assumption that the intrinsic geometry on E is regular even in the limit 6 - 
also assumed that the extrinsic curvature remains finite. 

Therefore, we obtain the following form of the Einstein-Hilbert action for the system including the singular hyper- 
surface E. 
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(8) 



where g± M „ is the induced metric, q± is the determinant of q± i iv, K± — q± K± fil/ is the trace of the extrinsic curvature 
K±^ v and <7±" is the inverse of q±^v In the expression (g) we have distinguished geometrical quantities in .M+ and 
M— by introducing the subscript ±, and have allowed the cosmological constant to have different values in these two 
regions. We have introduced the Lagrange multiplier field \^ v (y) to ensure the regularity of the intrinsic geometry of 
E. When we regularized the system and decomposed I eh into Iq and I± as in ([|), we implicitly assumed that the 
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induced metric and the extrinsic curvature are continuous across the boundaries B±. After taking the limit 8 — > +0, 
the extrinsic curvature remains finite but can be discontinuous across E. On the other hand, the induced metric should 
remain continuous across E even after taking the limit 6 — > +0 because of the finiteness of the extrinsic curvature. 
Provided that the hypersurface E is specified as the boundary of A4± by the parametric equation 

xf = Z£V), (9) 
the induced metric and the extrinsic curvature are given by 

q±nv(y) = e ±ti{y)e±Ay)g±MN\ x±=z±(y) , 

K±VLu{y) = g e ±u.{v) e ±Av)£n±9±MN\ x±=z±[y) , (10) 
where e± v are vectors tangent to E defined by 

r)7 M 

*Uv) = (ID 

and n± is the unit normal to E directed from A4_ to A4 + . To be precise, n+ is the inward-directed unit normal to 
E as the boundary of and n^ 1 is the outward-directed unit normal to S as the boundary of M_. 

Finally, the total action of the system is given by (Q), where I eh arid I ma tter are given by (|8|) and (||), respectively. 

III. VARIATIONAL PRINCIPLE 

In this section we derive equations of motion from the variational principle based on the action obtained in the 
previous section. Namely, we shall extremize the action I to t with respect to the variation 

9±mn{x) — > g±MN(x) + Sg±MN{x), 

Z^{y)^Z^{y) + 6Z^{y). (12) 

In the following we omit the subscript ± unless there is possibility of confusion. 

First, it is easy to show that the integrand of the volume term in Ieh changes as follows. 

y/fo\(R - 2A) - vTflf [(R - 2A) - (G MN + Ag MN )6g M N + (6g M % - 8g' M )., M + 0(S 2 )] , (13) 

where the semicolon represents the covariant derivative compatible with the background metric gMN (not with the 
perturbed metric gMN + SgMN ) , the indices M, N, ■ ■ ■ are lowered and raised by the background metric gMN and its 
inverse g MN , and Sg is defined by Sg = Sgjfi. Hence, 

MN i a „MN\ 



> I d D x ± y/\g~\(R - 2A) = - / d D x ±y /\g~\(G 

JM± J M± 



+ kg MN )5g MN 



Te I d^y^n™ {8g M . <N -5g., M )\ x±=z± 



Te f d D - l y^\n M 5Z¥ (R - 2A)| X±=Z± (14) 



The second term in the right hand side came from the total derivative in (13) and the last term is due to the change 
of the region to be integrated over. 

Next, let us consider the surface term in Ieh- As shown in ref. |p2f the variations of the induced metric and the 
extrinsic curvature are given by 

Sq^v = effe^(SgMN + SZm-,n + SZn-,m), 



§K v = ^n M n N (8g M N + 25Z m -n)K^ v 
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4 



where the right hand side is evaluated at x± = Z± (y) and 



(16) 



In order to make the covariant derivatives of SZ M well-defined, we have to extend 5Z M off £. The expressions ( |l5| ) 
are independent of the method of the extension. For details, see ref. p2|. Hence, 



q\K^^\q\K + 5(^\q\K)+0{S% 



where 



1 



-n L q^e*fe? (ST LMN + 5Z L , MN + R l , M lnSZ l '). 



Combining this with the second term in (flj), we obtain 



rN 



n M (SgZ-N ~ S9;m) + 25(y/\q\K)/y/\q\ = -(K^ - Kq^Sq^ - 2n M R MN 8Z 

' f /\q~\q^n M e?(6gMN + 25Z M ;N) 



where we have used the equations 



K,e„] M = 0, 
e "( e f n iv);M = eff(n N n N ). M = 0, 
q^ef?e» + en M n N =g MN . 

Thus, the variation of the Einstein-Hilbert action Ieh is given by 
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Now let us consider the variation of I m atter- 
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where q^ v is either q +fJ , u or <?_ p „, and 
V\gT\T^ N (x ± )=2 
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Therefore, SItot = is equivalent to the following set of equations. 
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r<MN , a MN ,.2r r MN 
— 9-jLlI/ = 0, 

Ki w _ Ri w = _ K 2 f s „ v _ _}_ Sq ^ j , _, 4 j 

^±iv = T< T ±MJV | x±=z± , (25) 

and 

a"" = -VW\(K V - K ^ v )- ( 26 ) 

In the right hand side of the last equation, the subscript — (or +) should be taken when Cq is written in terms of 
q +liv (or q-n„, respectively). 

Note that the equations (]24|) are the Einstein equation and Israel's junction condition pr|. The last equation is just 



to determine the Lagrange multiplier field X^ . Although the equation (25) looks like a new independent equation, 
it will be shown in the next section for simple examples that the equation is compatible with equations of motion of 
matter fields. Therefore, the action principle gives the correct set of equations: the Einstein equation, Israel's junction 
condition and equations of motion of matter fields. 



IV. SIMPLE EXAMPLES 



In this section we show that for simple examples, the equation (|25|) is compatible with equations of motion of matter 
fields. The first trivial example is the case in which all matter fields are confined on the hypersurface £. This case 
includes a shell with an arbitrary equation of state in a vacuum and the brane world scenario in a purely gravitational 
bulk with a bulk cosmological constant and arbitrary matter fields on the brane. In this case, the consistency condition 
( p5| ) is trivially satisfied since C± = and Cq does not change when Z± is changed with q^ fixed. 

As the second example, let us consider a simple case in which there is only a scalar field other than those matter 
fields confined on the hypersurface E. Namely, let us consider the following Lagrangian densities. 
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MN d M 3>d N <5> 



Co = £o(<fi+) + X<t>{^>+ - </>-), 
where Cq is the Lagrangian density for matter confined on S, and cf>± is the pullback of $ on £ defined by 

<t>±(y) = ®\x ± =Z ± {y) ■ 



(27) 



(28) 



The matter Lagrangian density Cq on E can depend on <f>+ as well. Note that the Lagrange multiplier field A^(y) is 
necessary in order that the scalar field should have single value on E and that E should be regular. For this example 
we can easily calculate T^' /JV , S^ v and F±m as follows. 
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where the right hand sides of 
Xa, — - 



j) and (31) are evaluated at x± — Z± (y), respectively. Hence, by using 



q\n 



;_=Z_ 



(29) 

(30) 
(31) 

(32) 



which is a part of equations of motion, it is confirmed that (Pq ) is satisfied. Thus, the consistency condition 
actually compatible with equations of motion of the scalar field. 

It is easy to extend the above analysis to an arbitrary number of scalar fields. 
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V. SUMMARY AND DISCUSSION 



We have presented an action principle of singular hypersurfaces in general relativity in any dimension without 



assuming any symmetry. Since an arbitrary number of scalar fields can be consistently included as shown in Sec. IV, 
the action principle is applicable to a wide class of scalar-tensor type theories of gravity in the Einstein frame. Besides 
the scalar fields, any kind of matter Lagrangian density on the hypersurface, which may depend also on the pullback 
of the scalar fields, can be included. The action principle is manifestly doubly covariant in the sense that coordinate 
systems on and off a hypersurface are disentangled and can be independently specified. More precisely, there are three 
independent coordinate systems: that on the hypersurface, those in two regions separated by the hypersurface. We 
have shown that, including variation of the metric, the position of the hypersurface and matter fields, the variational 
principle gives the correct set of equations of motion: the Einstein equation off the hypersurface, Israel's junction 
condition in a doubly covariant form and equations of motion of matter fields including the scalar fields. It is worth 
while mentioning that the position of the hypersurface measured from one side of the hypersurface and that measured 
from another side can be independently variated as required by the double covariance. 

Now let us discuss about application of the doubly covariant action principle to the brane world scenario. 

In refs. p4|-|30|] it was shown that the standard cosmology can be realized in the Randall-Sundrum brane world 
scenario in low energy as far as a spatially homogeneous and isotropic brane is concerned. After that, many authors 
investigated cosmological perturbations in the brane- world scenario |3l],^2|,[32] |4lJ . 

In particular, four independent equations for scalar perturbations on the brane in the plane symmetric (K = 0) 
background were derived recently by the author in rcf. (32). The number of independent equations is the same as 
in the standard cosmology, and it was shown that in low energy these sets of equations differ only by the non-local 
effects due to gravitational waves in the bulk. 

In the derivation of the four equations in ref. |32| the author took advantages of the doubly gauge invariant formalism 
developed in refs. [|T]j2^]. It was essential that the formalism includes perturbation of the position of a brane as a 
dynamical variable. Actually, as already discussed in ref. pH , if we fix the position of the brane by hand as in the 
Gaussian normal coordinate system, then it is in general inconsistent with convenient gauge choices in the bulk like 
a generalized Regge- Wheeler gauge []. In other words, as done in refs. p3] , |3"l| , we can construct D-gauge invariant 
variables from the perturbation of the position of the brane, and they are physical degrees of freedom independent 
of .D-gauge invariant variables in the bulk. The former gauge invariant variables are concise in the sense that it is 
localized on the brane, and the later variables can be expressed most concisely by the master variables introduced 
in ref. [ CT . Hence, the inclusion of the brane position as a dynamical variable provides us with the most concise 
configuration space. 

Now let us illustrate the above arguments about D-gauge- invariant variables by using some equations. For simplicity 
we consider perturbations around a background with 3-dimensional plane symmetry in 5-dimension. Namely, following 
the notation in ref. |p2j, we consider the metric 



and the imbedding relation 



ds\ = g M Ndx M dx N = (gfj N + 5g M N)dx M dx N (33) 



x M = Z M {y) = Z^ M {y) + 5Z M (y), (34) 



where the background is specified by a plane-symmetric background metric 

3 



g ( M N dx M dx N = lab dx a dx b + r 2 Y^{dx 1 ) 2 (35) 



i=i 



and such background imbedding functions Z^ M (y) that Z^ a depend only on y° and that Z^ 1 = y l . Here, the 
two-dimensional metric j a b and the function r 2 are assumed to depend only on the two dimensional coordinates {x a }. 
As for perturbations, since in the linear order the perturbations of the position of the hypersurface are decoupled 
from vector and tensor perturbations, we consider scalar perturbations: 



2 In the literature it is sometimes called a generalized longitudinal gauge. 
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Sg M Ndx M dx N = J d 3 k [h ab Ydx a dx b + 2h (L)a V (L)t dx a dx l 
+(h(LL)T(LL)ij + h( Y )T{Y)ij)dx l dx ] ] , 
5Z M dx M = J d 3 k [z a Ydx a + z^Vwidx*] , (36) 

where Y = cxp(— ik ■ x), V{L)i — d{Y, T(LL)ij = 2<9 i d 7 F + (2k 2 /3)<5y-Y and T(Y)ij = SijY, and all coefficients are 
supposed to depend only on the 2-dimensional coordinates {x a } of the orbit space. Here, x denotes coordinates {x 1 } 
of the three-dimensional plane (i = 1, 2, 3), and k represents the momentum {fc,} along the plane. Hereafter, we omit 
k in most cases. It is easy to see how the coefficients {h's, z's} transform under the 5-gauge transformation and to 
construct 5-gauge-invariant variables. Therefore, we obtain the following 5-gauge-invariant variables. 

4>a = z a + X a , (37) 

and 

Fab = h a b — V a Xi, — VbX a , 

2k 2 

F = h (Y) - X a d b r 2 + h (LL) , (38) 

where X a = htL)a ~ T 2 d a {r~ 2 h(j, r\) and V a represents the covariant derivative compatible with the 2-dimensional 
metric 7 a &. The former variables ( p7| ) correspond to perturbations of physical position of the hypersurface S and its 
normal component (j) a n^ a appears in the doubly-gauge-invariant junction condition, where 

n (0)a ig 

the background 

unit normal to the hypersurface. The latter (p8|) correspond to gravitational perturbations in the bulk and can be 
most concisely expressed in terms of the master variable $ as 

Fab = ~ r (V a V 6 $ - ^V 2 $ 7ab + 7^*7afc 

F= T - (V$ - . (39) 

The perturbed Einstein equation in the bulk is reduced to the following simple equation called master equation: 

r 2 V a [r _1 V a (r~ 1 $)] - k 2 r~ 2 $ = 0. (40) 

In the generalized Regge- Wheeler gauge where h^ a = ^(ll) = 0, the 5-gauge-invariant variables are given by 
4>a = z a , F ao — h a i, and F = hry)- On the other hand, in the Gaussian normal gauge where z a = zij\ = h(L)a n — 
h a bn^ b — 0, these are given by cf> a — X a and (|38|). Note that in the Gaussian normal gauge, <j> a is expressed in terms 
of metric perturbation. Therefore, it is evident that <p a cannot be set zero even in the Gaussian normal gauge since 
X a in general. Actually, requiring <f> a = in the Gaussian normal gauge is equivalent to requiring z a — in the 
generalized Regge- Wheeler gauge, which is not possible in general. 

Of course, it is always possible to take the Gaussian normal coordinate system. In this coordinate system, as 
illustrated above, the 5-gauge-invariant variable <f> a is expressed in terms of metric perturbations. Hence, as done in 
ref. |Q for a static background by a gauge-dependent method, we need to extract degrees of freedom of a n^°^ a from 
the metric perturbations. Classically, this procedure should not be difficult since we can use equations of motion. 
However, quantum mechanically, we have to be careful when we use the equations of motion to reduce the action. 

The next task in the future is to obtain the second-order action for perturbations by using the doubly covariant 
action obtained in this paper. After that, we need to obtain the corresponding reduced action by using a formalism to 
treat constrained systems, eg. Dirac's method |^3| or Faddeev-Jackiw method Q]. As shown in ref. [Q, perturbative 
behavior of the Wheeler-de Witt wave function can be investigated by the usual quantum field theory in curved 
spacetime with the reduced action. 
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